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ABSTRACT 

The nature of the progenitor systems of Type la supernovae is still uncertain. 
One way to distinguish between the single-degenerate scenario (SDS) and double- 
degenerate scenario (DDS) is to search for the post-impact remnant star. To 
examine the characteristics of the post-impact remnant star, we have carried 
out three-dimensional hydrodynamic simulations of supernova impacts on main 
sequence-like stars. We explore the evolution of the post-impact remnants using 
the stellar evolution code MESA. We find that the luminosity and radius of the 
remnant star dramatically increase just after the impact. After the explosion, 
post-impact companions continue to expand on a progenitor-dependent timescale 
of ~ 10 2 ' 5-3 yr before contracting. It is found that the time evolution of the 
remnant star is dependent not only on the amount of energy absorbed, but also on 
the depth of the energy deposition. We examine the viability of the candidate star 
Tycho G as the possible remnant companion in Tycho's supernova by comparing 
it to the evolved post-impact remnant stars in our simulations. The closest 
model in our simulations has a similar effective temperature, but the luminosity 
and radius are twice as large. By examining the angular momentum distribution 
in our simulations, we find that the surface rotational speed could drop to ~ 
10 km s _1 if the specific angular momentum is conserved during the post-impact 
evolution, implying that Tycho G cannot be completely ruled out because of its 
low surface rotation speed. 

Subject headings: binaries: close, — stars: evolution, — methods: numerical, - 
supernovae: general, — supernovae: individual (Tycho's SN) 
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INTRODUCTION 



Type la supernovae (SNe la) are important as cosmological distance probes and 
sources of metal enrichment for the interstellar medium. Although they can b e treat ed 



as "s t andardizable candles," intrins ic variations in SN la properties do exist (ILivio 



2000 



Hillebrandt &: Niemeyer 



2000l ). In principle these could be associated with 



different SN la progenitor systems and explosion mechanisms. It is widely accepted 



that SNe la are thermonuclear explosions of carbon-oxygen (CO) white d warfs 



fWDs) 



but so far observations have not conclusively favored only a single mo del (jLivio 



Hillebrandt fc Niemeyer 



2000 



Ruiz-Lapuente 



2012 



Wang fc Han 



20121 ) 



2000 



Whelan fc lb en 



The mainst r eam p r ogenitor 



Iben &: Tutukovj ( 



1973 



198J) 



Nomoto 



scenar ios include the single-degenerate scenario (SDS; 



Webbink 



19821)) and the double-degenerate scenario (DDS; 



1984 )). The SDS involves a CO WD accreting mass 



from a non-degenerate companion, such as a main sequence (MS) star, red giant (RO), 
or helium (He) star, through Roche-lobe overflow (RLOF). Once the CO WD reaches the 
Chandrasekhar mass, it becomes unstable and explodes as an SN la. The DDS instead 
involves gravitational wave-driven merging of two CO WDs with total mass greater than 



the Chandrasekhar mass. Recent ob servations of the nearby SN 2011 f e provide use 



constraints on its progenitor system 



Brown et al. 



2012 



Chomiuk et al. 



Li et al. 



2011 



Nugent et al. 



2011 



Bloom et al. 



ul 



2012 



20121 ) . though both scenarios may contribute to the 



SN la population. 

To distinguish between the SDS and DDS, different approaches have been proposed 
during the past decade. For example, detection of the pre-supernova circumstellar medium 
around supernova remnants (SNRs) provides evidence of the mass transfer phas e in the SDS 



( iBranch et al. 



1995 



Sternberg et al. 



2011 



Foley et al. 



2012 



Dilday et al. 



2012|). Collision 



of the SN ejecta with a non-degenerate companion may affect the early SN light curve, 
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depending on the companion type and viewing angle (lKasenll2010i ). Numerical simulations 
of the SDS also suggest that the non-degenerate companion will survive th e SN impact anc 



shou 



2012|). 



d be detectable in the SNR ( Marietta et al. 



2000 



Pakmor et al. 



2008 



Pan et al 



2010L 



Ruiz-Lapuente et al. 



(120041 ) found a star, which they called Tycho G, that is located 
in the central region of Tycho's SN and has a high proper motion, which in the SDS 
would be associated with the original orbital s peed. This sugges t s that Tycho G could 



Kerzendorf et al. 



(120091 ) found that the 



be a post-impact SDS remnant star. However, 
upper limit of the rotational speed of Tycho G is < 7.5 km s" 1 , which would seem to be a 
problem for the SDS inte rpretation, since the companion stars in close binary systems are 



usually rapidly rotating. 



Pan et al. 



( 120121 ) address this problem by analyzing the angular 



momentum lost during the supernova impact, but the final rotational speed depends on the 



binary evolution on non-degenerate companions have been 1 


Dublished ( 


Hachisu et al. 


1999: 


Han & Podsiadlowsk 


1 2004 


Han 


2008; 


Hachisu et al. 2008b 


Wane: et al. 


2009 




Mens: et al. 


2009; 


Meng & Yang 


2010 




Wang & Hani 


2010 


), usually they have not included the SN 



impact and post-impact evolution. 



Podsiadlowskil ( 120031 ) studied the post-impact evolution of a 1M subgiant by assuming 
a fixed amount of energy input and mass stripping. Unsurprisingly, they found that the 
post-impact star could be overluminous or underluminous compared to a normal star of 
the same tempera t ure, d epending on the choices of these parameters. A more recent study 



by 



Shappee et al. 



(120121 ) examined the post-impa ct evolution of a 1M MS star from the 



hydrogen cataclysmic varia bles (HCV) scenario i n 



Marietta et al. 



(120001 ). using an approach 



similar to that adopted by iPodsiadlowskil (120031 ) . They find an overluminous post-impact 
companion (10 — 1O 3 L ) and suggest that because of this Tycho G may not be associated 
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with Tycho's SN. However, without detailed study of the impact of the SN ejecta on the 
companion star, especially shock compression in the stellar interior and the depth of energy 
deposition, these results cannot be considered definitive. Furthermore, the stellar mass and 
delay time for Tycho G might differ from the p rogenitor in the HCV scenario, which is the 



only model considered in 



Shappee et al. 



tenia ) 



In this paper, we consider a wide 
channel (WD+MS or WD+subgiant) in 



range of progenitor m odels from the supersoft 



Hachisu et al. 



(l2008bl ) (hereafter HKN) and perform 



three-dimensional hy drodynami c s sim ulations of the impact of SN la ejecta using the 



method described in 



Pan et al. 



(120121 ) (hereafter PRT). We then use a one-dimensional 
stellar evolution code to study the post-impact evolution for each of these progenitor 
models. In the next section, the numerical methods and progenitor models are described. 
In Section EJ we present the results of post-impact evolution for different companion 
models and compare their observational quantities with Tycho G. In the final section, we 



summarize our results and conclude. 



2. NUMERICAL METHODS 

Our numerical simulations are divided into three stages. The first stage uses a 
one-dimensional stellar evolution code to construct the progenitor system at the onset of 
the SN la based on detailed binary evolution. The second stage uses a three-dimensional 
hydrodynamics code to simulate the impact of SN la ejecta on the binary companion. The 
final stage takes the post-impact companion star from the hydrodynamics simulation as 
input to the stellar evolution code, which is then used to simulate the post-impact evolution. 
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The hydrodynamics code used is FLAS 



2.1. Simulation Codes 



version 3 (IFryxell et al 



2000 



Dubey et al. 



20081 ). FLASH is a parallel, multi-dimensional hydrodynamics code based on block- 



structured adapt ive mesh refinement (AM R). The equation of state (EOS) applied is the 



Helmholtz EOS ( Timmes fc Swesty 



20001 ) . The Poisson eq uation is soly ed on the AMR 



20081 ). The detailed 



mesh using a Fourier transform-based multigrid algorithm (IRicker 
numerical setup and initial conditions are the same as in PRT. In this paper, we consider 
only three-dimensional simulations that include orbital motion. The i nitial binary systems 



are assumed to be in RLOF, and the SN model used is the W7 model (INomoto et al 



1984h . 



To construct the progenitor systems and simulate the post-impact evolution of the 
remnant star, the stellar evolution code MESA^ (Modules for Experiments in Stellar 



Astrophysics; 



Paxton et al. 



(120111 )) is e mployed. The mixin g-length theory (MLT) of 



convection used is the modified MLT by 



Henyey et al. 



(119651 ). which allows the convective 



efficiency to vary with the opacity and has important effects on convective zones near the 
surfaces of stars. The ratio of mixing length to the local pressure scale height is set to 
a = 1.918, which is the value calibrated using the Sun. The initial metallicity is assumed 
to be Z = 0.02 for all the progenitor models. 



2.2. Progenitor Systems 

The progenitor systems are taken from the models in HKN. HKN studied the binary 
evolution of WD+MS systems and found the region of the donor mass-orbital period 
plane for which SNe la may occur. In addition to RLOF mass transfer, HKN consider the 

1 http : / /flash . uchicago . edu 
2 http : //mesa. sourcef orge .net 
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optically thick wind from the WD and the mass stripping effects of a massive circumstellar 
torus. Although rece nt EVLA observations of SN 20 life disagree with this scenario 



( jChomiuk et al 



20121 ) . HKN can explai n the circumste l 



SN 2005gj, and SN 2006X. Furthermore, 



Hachisu et al. 



ar matter around SN 2002ic, 



(j2008al ) found that the theoretical 



delay-time distribution (DTD) based on this scenario is consistent with the observed DTD. 

Figure 7 and Figure 8 in HKN show the initial and final regions of the donor 
mass-orbital period plane for WD+MS systems with initial WD mass Mwd = 1.0M Q , 
hydrogen composition X = 0.70, metallicity Z = 0.02, and stripping rate parameter 
Ci = 3. We take six evolved MS models from these figures as our progenitor models (see 
Table [1]). Since these models are assumed to be in RLOF, we use MESA to evolve a 
zero-age-main-sequence (ZAMS) star with the initial mass and metallicity taken from HKN 
until the radius of the MS star reaches its Roche-lobe radius. The Roche-lobe radius, Rl, 
is approximately given by 



Eggletonl fl!983h : 



Rl 

a 



0.49g 2/3 



1) 



0.6g 2 / 3 + ln(l + gV3) ' 

where q is the mass ratio, and a is the binary separation. The binary separation can be 
found using Kepler's third law with the given mass and period from HKN. Although MESA 
has the ability to model binary evolution, we do not include the detailed physics of binary 
evolution as HKN considered but rather assume the MS star to have a constant mass-loss 
rate. We vary the mass-loss rate until the MS star matches the final mass and radius given 
in HKN within a timescale comparable to theirs (~ 10 6 yr). The mass and radius of these 
MS stars are shown in Figure [TJ Since we assume RLOF, the ones with longer periods have 
larger radii and longer delay times. 

Table [T] summarizes the initial and final conditions of these six models we have created 
using MESA. The "MS-r" case in PRT is also included as a c omparison (the star G in 



Table [Q. Note that the HCV model in 



Marietta et al 



fl200Q[ ): 



Shappee et al. 



(120121 ) and 
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the HCV model in PRT are not considered here, since these scenarios did not consider 
detailed binary evolution. Figure |5] shows the density profiles of all these models at the time 
of the explosion. Most models are slightly evolved MS stars, but some of them are close to 
ZAMS, depending on the initial periods in HKN. MESA provides all the thermodynamical 
quantities we need for FLASH simulations. To save computation time, the compositions of 
the companion models are adjusted to 1 H, 4 He, plus 12 C, where 12 C represents all elements 
heavier than helium. After converting MESA models to FLASH models, we perform 
FLASH simulations using the method described in PRT. 

2.3. Post-Impact Companion Models 

After the SN impact, the companion star is heated and loses about ~ 10 — 20% of 
its mass due to stripping and ablation by the SN ejecta (see PRT for detailed description 
of the SN impact). The remaining mass and final conditions of the remnant stars are 
summarized in Table |2j Figure [3] shows a typical gas density distribution for model F in 
the orbital plane at the end of the simulation (t = 3.26 x 10 4 sec). It is difficult to run 
hydrodynamics simulations up to the age of historical SNRs (~ 500 — 1000 yrs) because 
the dynamical time scale, which determines the timestep, is only on the order of hundreds 
of seconds. Moreover, the star is close to being spherically symmetric. Therefore, the most 
straightforward method is to turn the post-impact results back into a one-dimensional 
problem and solve for the subsequent evolution using MESA. 

To run simulations of the post-impact evolution of remnant stars in MESA, the 
first step is to convert the three-dimensional data into angle-averaged one-dimensional 
radial profiles. Our setup in MESA cannot deal with the spin of the remnant star, and 
the angular momentum is finite at the end of the FLASH simulations. However, PRT 
demonstrated that the loss of angular momentum during the SN impact will significantly 
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Table 1: The progenitor models 



Model 


M a 


Pi 


Ri 






M f b 


Pi 


Rt 


Lt 


T e S,{ 




(M e ) 


(day) 


(Re) 


(L Q ) 


(K) 


(M Q ) 


(day) 


(Re) 


(L Q ) 


(K) 


A 


2.51 


0.477 


1.83 


39.2 


10,696 


1.88 


0.350 


1.25 


2.35 


6,392 


B 


2.51 


0.600 


2.08 


42.4 


10,224 


1.92 


0.466 


1.50 


3.64 


6,516 


C 


3.01 


1.23 


3.64 


110.0 


9,800 


1.82 


1.09 


2.63 


8.06 


6,003 


D 


2.09 


0.472 


1.67 


19.2 


9,358 


1.63 


0.353 


1.19 


2.09 


6,372 


E 


2.09 


0.589 


1.91 


20.8 


8,933 


1.59 


0.470 


1.42 


3.15 


6,450 


F 


2.09 


0.936 


2.59 


23.9 


7,934 


1.55 


0.770 


1.97 


5.09 


6,182 


G* 


2.00 


1.00 


1.70 


17.6 


9,083 


1.17 


0.233 


0.792 


0.463 


5,355 



a The mass (Mi), period (day), radius (Ri), luminosity (Li), and effective temperature (T e ff,i) for different 
progenitor models at the beginning of RLOF for WD+MS systems, using the initial masses and orbital 
periods in Figure 7 of HKN. 

& The mass (Mf), period (day), radius (Rf), luminosity (if), and effective temperature (T e gj) for different 
progenitor models at the time of the SN explosion, using the final masses and orbital periods in Figure 8 of 
HKN. 

*The MS-r model in PRT. 
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Fig. 1. — Stellar mass and radius for all considered progenitor models at different stages: ini- 
tial state (x), just prior to SN (filled circles), after reaching hydrostatic equilibrium following 
SN (filled squares). 
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10 3 




m/M, 

Fig. 2. — Density vs. radius profiles for all the MS models at the time of the SN explosion. 
The label for each model refers to Table [D 
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Table 2: Post-impact remnant stars 



Model M SN (M ) a i? SN (R Q ) AM/M, w Linear (km s" 1 ) b L SN (£©) T cff , S N (K) 



A 


1.64 


3.87 


12.8 


% 


179.0 


16.9 


5,954 


B 


1.65 


4.76 


14.1 


% 


178.7 


22.4 


5,760 


C 


1.56 


7.61 


14.3 


% 


136.3 


44.2 


5,398 


D 


1.43 


3.42 


12.3 


% 


187.7 


13.1 


5,936 


E 


1.44 


3.91 


9.4 


% 


190.9 


14.6 


5,715 


F 


1.30 


4.09 


16.1 


% 


142.6 


15.0 


5,618 


G 


0.93 


4.45 


20.1 


% 


271.0 


13.9 


5,289 



a The mass (Msn), radius (Rsn), percentage of unbound mass (AM/M*), linear velocity (^Linear), luminosity 
(Lsn)i arid effective temperature (T e ff,SN) of initial post-impact hydrostatic models in MESA. 
^Linear velocity includes the pre-supernova orbital speed and kick velocity. 
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Fig. 3. — Gas density distribution in the orbital plane for a three-dimensional SN la simu- 
lation with model F (see Table H]) at time t = 3.26 x 10 4 sec after the explosion. The frame 
shows a domain spanning 10R & . The color scale indicates the logarithm of the gas density. 
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decrease the spin. Since the rotational velocity is significantly less than breakup, 
Qk = (GM/Rl) 1 / 2 ~ 5 x 10~ 5 s _1 , we ignore the spin during MESA simulations, but 
consider it separately using the post-impact specific angular momentum from FLASH. The 
role of rotation will be discussed in more detail in Sectior j373l 

Since the companion models at the end of the FLASH simulations are not yet in 
hydrostatic equilibrium, the averaged radial profiles of density, temperature and pressure 
cannot be used directly in MESA. However, the averaged specific entropy and composition 
mass profiles are conserved if the system is adiabatic. Since the Mach numbers at the end of 
the FLASH simulations are subsonic for most gravitationally bound gas, we can use these 
profiles, as a start, to construct hydrostatic models for MESA. Hydrostatic equilibrium 
should be achieved relatively quickly after the SN impact in comparison with the thermal 
timescale. 

To calculate the averaged specific entropy and composition profiles, only gravitationally 
bound gas is considered, and the center of the remnant star is taken to be the location of 
the gravitational potential minimum. The simulation domain is then divided into 128 — 256 
radial bins, depending on the zone spacing and the distance between the star's center and 
the boundary of the simulation box at the end of the FLASH simulation. Although FLASH 
and MESA include a more realistic EOS, for simplicity only an "ideal gas plus radiation" 
EOS is used to calculate the entropy profile. Therefore the specific entropy is 



where Na is the Avogadro constant, is the Boltzmann constant, a is the radiation 
constant, p is the bound gas density, T is the temperature, and /x is the mean molecular 
weight. Using the total bound mass (M*), the EOS, and the entropy and composition 
profiles, the density, pressure (P), temperature (T), and radius can be calculated by solving 




(2) 
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the hydrostatic equations, 

(3) 



dM r 47rr 2 p' 
dP GM r 



(4) 



dM r 4vrr 4 ' 

where G is the gravitational constant, and M r is the mass within radius r. The hydrostatic 
equations are solved using the fourth-order Runge-Kutta method with adaptive stepsize 
control. We use the shooting method, varying the initial central density to match the 
boundary condition M(r = R*) = M*, p(r = R*) = 0, and P(r = R*) = 0. 

The reconstructed mass and radius of the post-impact remnant stars are shown in 
Figured! The upper left and right panels in Figure H] and Figure show the averaged ID 
entropy and helium composition profiles from the FLASH output. Models B, C, E, and 
F have larger radii and higher central helium abundances (upper right panel), since their 
delay time and initial orbital periods are longer in HKN. We flatten the entropy profiles in 
the outermost region (0.95 — 0.99 < m/M* < 1) to avoid negative entropy gradients. The 
sensitivity of our results to these flattened entropy profiles have been tested by varying the 
amount of the flattened entropy and using a smoothly increasing profile instead of a flat 
one. We find that the outermost region of the entropy profile does not significantly affect 
the hydrostatic solution and post-impact evolution. 

Another necessary variable is the initial luminosity profile, L(m). If there is no 
convection, the luminosity profile can be estimated using the radiative temperature gradient 
expression, 

L (™) = o > (5) 

6k dm 

where k is the opacity and c is the speed of light. The opacity can be calculated using 
the kap module in MESA for a given density, temperature, and composition. This 
estimate is reasonable only at the beginning, since the modified entropy gradient is positive 
everywhere. The energy deposited by SN ejecta heating causes a small temperature bump 
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at m ~ 0.95M* (see middle panels in Figure H] and Figure [5]), leading to the luminous region 
in the luminosity profiles in Figure H] and Figure EJ 

3. EVOLUTION OF THE POST-IMPACT REMNANT STAR 

In this section, the post-impact evolution of remnant stars using MESA is qualitatively 
described and compared with Tycho G. The effect of post-impact evolution with a different 
explosion energy is also studied to provide an understanding of the resultant evolution 
for a non-W7 explosion. In addition, the angular momentum distribution of post-impact 
remnant stars is examined by using the specific angular momentum from FLASH and the 
hydrostatic models in MESA, addressing the rotation problem in Tycho G. 

3.1. Post-Impact Evolution 

The hydrostatic solutions described in the previous section are used as initial conditions 
for the models used in MESA. The initial timestep in MESA is chosen to be between 
10 -4 and 10" 1 years, depending on the particular model. The compositions in the 
FLASH simulations are expanded to the extended network in MESA with 25 isotopes. 
The percentages of compositions other than X H and 4 He are scaled to solar abundances, 
but with the same total metallicity. When starting a MESA simulation the model is 
iterated until convergence to ensure that both the EOS and the energy generation rate are 
consistent. During the first few steps, the evolution is transient, but eventually evolves to a 
parameter-independent model within a year. Thus, only simulation results after one year 
are considered in this study. 

Figure [6] shows the evolution of the photospheric radius, luminosity, and effective 
temperature as functions of time. The post-impact remnant stars rapidly expand on a 
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Fig. 4. — Hydrostatic solutions for entropy (S) helium composition (Y), density (p), tem- 
perature (T), pressure (P), radius (R), and enclosed luminosity L(m) for models A, B, C, 
and D in Table [2j 
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Fig. 5- 



Similar to Figure H] but for models E, F, and G. 
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Table 3: Post-impact remnant star after 440 years 



Model 


M SN (M ) 


R SN (Rq) 


£sN (L & ) 


^eff,SN (K) 


^Rotation (km S 


A 


1.64 


7.18 


202.7 


8,135 


10.5 ± 1.3 


B 


1.65 


11.9 


279.0 


6,832 


8.3 ±0.6 


C 


1.56 


10.5 


81.4 


5,356 


2.0 ±2.9 


D 


1.43 


6.57 


111.4 


7,319 


9.1 ± 1.3 


E 


1.44 


4.56 


20.2 


5,737 


26.7 ±4.4 


F 


1.30 


4.35 


17.0 


5,623 


17.3 ±2.2 


G 


0.93 


4.37 


13.4 


5,287 


18.9 ±2.5 


GH09 t 


1.0 


1.4 - 2.8 


1.9 - 7.6 


5900 ± 100 




GH09 t 


1.4 


1.6 - 3.3 


3.0 - 11.8 


5900 ± 100 





a The mass (Msn), radius (Rsn), luminosity (Lsn), effective temperature (Ttf^sN), and surface rotational 
speed (^Rotation) of post-impact remnant stars at 440 yr after SN explosion. 

T The radius and luminos ity are estimated using the surface gravity value log g = 3.85 ± 0.3 in 



Gonzalez Hernandez et al 



(|2009i) . assuming the mass is 1M© or 1.4M© 
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timescale of ~ 10 2 years for models A, B, D, and E, and ~ 10 3 years for models C and F. 
Note that there is no significant change in model G within 10 4 years, because the thermal 
timescales in model G are much longer than other stars. 

Here we describe the detailed evolution of model A as an example. Due to the high 
opacity in the outermost region of the envelope, the photospheric luminosity is ~ 4 L , but 
the luminosity inside the envelope is much higher than the surface luminosity (Figure HJ). 
Therefore, the strong radiation at m ~ 0.9M* expands the outermost < 1% of mass. In 
addition, the luminosity profile of the outermost 10% of mass flattens due to radiative 
diffusion. This radiative d iffusion timescale is characterized by the local thermal timescale 



( jHenyey fc L'Ecuyer 



19691 ). 

Tth(/ 



647TCT' 



SB 



J"3 



pdr 



(6) 



where age is the Stefan-Boltzmann constant, and C p is the specific heat capacity. Although 
the global thermal timescale, £ t h = GM 2 /2RL, is of the order of 10 6 years, the local thermal 
timescale in the envelope region is only ~ 10 2 — 10 3 years. 

After ~ 200 years, the deposited energy has been radiated away and model A starts 
to contract by releasing gravitational energy. The luminosity decreases, but the effective 
temperature increases. The star will eventually return to a stage similar to the ZAMS on 
a global thermal timescale. As an SN la remnant may not be recognizable in millions of 
years, the simulation is terminated at 10 4 years. 

The evolutionary tracks of post-impact remnant stars in the Hertzsprung- Russell (H-R) 
diagram are plotted in Figure [7J In general, the remnant stars become brighter but cooler 
after the impact. By comparing Figure [7J with Figure [TJ we find that there is a trend for 
more massive stars coupled with smaller radius (shorter orbital period) to evolve faster in 
the HR diagram. In addition, stars with smaller radii have higher effective temperatures 
after the SN impact. This can be understood given the fact that stars with shorter orbital 
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Fig. 6. — Evolution of the photosphere radius, luminosity, and effective temperature as 
functions of time. Each line shows the evolution of a post-impact companion star in Table |2j 
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periods are closer to the exploding WD, leading to more violent SN impacts and greater 
energy deposition, but with less unbound mass due to their more compact states. 



The Tycho G star, b ased on observations 



also plotted in Figure [0 



Ruiz-Lapuente et al. 



by 



Gonzalez Hernandez et al. 



()2009|), is 



(120041 ) determined Tycho G's effective 
temperature, T eff = 5750 ± 250 K, and surface gravity, log (q /cm s~ 2 ) = 3.5 ± 0.5, by 



fitting spectral lines. Subsequently, 



Gonzalez Hernandez et al. 



(120091 ) updated the observed 



values to T e ff = 5900 ± 100 K an d log(g/cm s 2 ) = 3.85 ± 0.30 based on fits to iron lines. 



Gonzalez Hernandez et al. 



( 120091 ) estimated the radius of Tycho G using the surface gravity, 
assuming the mass of Tycho G to be 1M . They estimated the bolometric luminosity 
to lie in the range 1.9 < L*/L Q < 7.6. However, this luminosity could be underestimated 
if Tycho G is more massive. Since log oc 2 log R + 4 log T and log g oc log M — 2 log R, 

could be increased to 3.0 < L*/L < 11.8 if the mass of Tycho G were instead 
1.4M . Our post-impact remnant stars have masses ranging from 1.3M to 1.65M . A 
mass M* = 1.4M instead of 1M is, therefore, assumed for Tycho G. As a result, the 
placement of Tycho G in the HR diagram is now closer to our models A, B, D, E and 
F immediately after the SN impact (triangle symbols in Figure [7]). However, Tycho's 
SN exploded 440 years ago, and models A, B and D would have evolved to a hotter and 
more luminous state in that time (star symbols in Figure Ej), suggesting that model E 
among our progenitor systems is the least discrepant with Tycho G. Our model E has a 
consistent effective te mperature, but the luminos i ty (ra dius) is twice brighter (larger) than 



the observed value in 



Gonzalez Hernandez et al. 



(120091 ) . Model F has a luminosity similar 



to that of model E, but the effective temperature is a few hundred degrees lower than 
Ty cho G. A less massive or smaller model than model E may better match the observation 



by 



Gonzalez Hernandez et al. 



(120091 ) . Table |3] summarizes the stellar properties for all the 
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progenitor models at 440 yr after the SN explosion! 



3.2. Dependence on the SN Ejecta Energy 



Although Tycho' s SN has previously been su ggested to be a sublu minous SN la 



(Ivan den Berg 



19931 ). or a overluminous SN la flRuiz-Lapuente 



200J), the most 



recent observation using th e light echo suggests that Tycho's SN is actually a normal 



SN la ( IKrause et al 



20081 ). In our hydrodynamics simulations, the W7 model 



(i?SN = 1-23 x 10 51 erg) is assumed to describe the SN explosion. However, the exact 
explosion energy could be different from case to case. Therefore, two additional simulations 
for models D and E, which are the remnant stars closest to Tycho G, are performed with 
double the explosion energy. We specify that this additional energy is in a thermal form. 



Pakmor et al. 



(120081 ) have shown that the unbound mass of the companion star is 
linearly proportional to the SN kinetic energy. In our hydrodynamics simulations, a similar 
trend is observed, but it is not linear, since there are essentially two mechanisms that 
unbind the mass: ablation and stripping (see PRT for detailed description). The resulting 
mass of the post-impact remnant star for model D (model E) is Mw7 = 1.43 M Q (1.44 M & ) 
for the W7 case and M 2 x = 1-23 M & (1.37 M ) for the case with double the explosion 



Two recen t obse rvations of SN 1006 by 



Gonzalez Hernandez et al. 



Kerzendorf et al. 



(120 12 ) and 



( 120121 ) suggest that there are no surviving evolved companions in the cen- 
tral region of SN 1006. We have noted that the star B90474 (\ogg = 3.05 ± 0.12 (cm s~ 2 ), 
T eff = 5051 ± 38 (K)) and B14707 (logs = 3.36 ± 0.15 (cm s" 2 ), T cfr = 5065 ± 47 (K)) have 
similar surface gravity as our Model G (logg = 3.14 (cm s -2 ), T g = 5288 (K)) but a lower 
effective temperature. However, the star B14707 lies much closer than the SNR and the 
distance uncertainity of star B90474 is large. 
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Fig. 7. — The evolutionary tracks in the HR diagram for different post-impact companion 
models. Each line shows an evolutionary track of a companion model in Table [2] over an 
interval of 10 4 years. The filled circles indicate the condition of stars just before the SN la 
explosion; filled triangles present the condition of stars at ~ 1 year after the SN impact; 
star symbols show the conditions at 440 years after the SN impact, which is the same age 
as Tycho's SN; filled squares present the conditions at 3, 000 years after the SN impact. 
The star symbols with er ror bars show the observed l umino sity and effective temperature of 



Tycho G as measured by 



Gonzalez Hernandez et al. 



(|2009|) (GH09) in Table El 
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energy. 

Figure [H] shows the temperature profiles of these four explosion cases after the SN 
impact. The heating is stronger in the higher explosion energy case, resulting in a lower 
remnant mass and causing a lower central temperature. The temperature bump in the 
higher explosion energy case shifts inward because the heating penetrates to deeper mass 
layers with a stronger impact. Therefore, the radiative diffusion time for the propagation of 
energy to the surface becomes longer. 

The evolution of the photospheric radius, luminosity, and effective temperature can 
be seen in Figure The lower mass and temperature in the higher energy case cause the 
radius, luminosity, and effective temperature to be smaller than for the W7 case. In general, 
increasing the explosion energy will lower the effective temperature and luminosity, but the 
evolution timescale is also proportional to the explosion energy. For example, a dramatic 
difference in the evolution is seen for model D in which the luminosity decreases by a factor 
of eight (see Figure [9]) at 440 years in comparison with the standard explosion energy (W7). 
This difference reflects the deeper energy deposition and longer thermal diffusion time 
associated with increased explosion energy. If we desire to match the luminosity of Tycho 
G by enhancing the explosion energy, the effective temperature becomes too low since the 
mass of the remnant star becomes lower. Thus, neither an overluminous nor a subluminous 
explosion model for model D or E can perfectly match the observed properties of Tycho G. 



3.3. Rotation Problem 



In the binary evolution scenarios in HKN, the evolution time from the onset of RLOF 
in the MS+WD systems to the SN la phase is about 10 6 years. Although the mass t ransfe r 



in close binary systems may cause the companion stars to be non-corotating, 



Zahnl fll977h 
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Fig. 8. — Temperature vs. mass profiles of the hydrostatic solutions for model D (thick lines) 
and model E (thin lines) with different SN la explosion energies. The solid lines indicate 
the temperature profiles with explosion energy using the W7 model (-EgN = 1-23 x fO 51 erg). 
The dashed lines show the temperature profiles with double the explosion energy. 
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Fig. 9. — Similar to Figure El but for model D (thick lines) and E (thin lines) with different 
SN la ejecta energy. 
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finds the synchronization time in binary systems to be 

4 y nc ~ q~\a/Rf ~ 10 4 ((1 + q)/2q) 2 P A years, (7) 

which is ~ 10 3 years in our cases. Therefore, the companion stars should all be synchronized 
by tidal locking and be rapidly rotating, resulting in surface rotational speeds of the order 
of hundreds of kil ometers per second. How ever, the upper limit on the rotation speed of 



Kerzendorf et al. 



(120091 ) is v sinz < 7.5 km s , where i is the inclination 



Tycho G found by 
angle. 

PRT indicate that the companion star (model G in Table [T]) loses about half of its 
angular momentum while only losing about 20% of its mass. With the strong assumptions 
that the parameter a is constant in the angular momentum expression, J = aMR 2 u, and 
that the companion star is in a state of solid-body rotation before and after the impact, 
the rotation speed can significantly drop to about a quarter of its original rotation speed. 
However, the equilibrium status of the remnant stars and the location of the photosphere 
were not calculated in detail in PRT. Combining the post-impact angular momentum 
distributions in FLASH with the stellar structure of post-impact remnant stars in MESA, 
the rotation problem in the Tycho G can be addressed quantitatively. 

Before the SN impact, the companion stars were set into uniform rotation and were 
given a spin with spin-to-orbit ratio of 0.95. After the SN impact, the companion stars 
are no longer in solid-body rotation. The angle-averaged angular velocity profiles can be 
calculated using 

n(r)-/ V4> 



. r cos 9 1 <f>,e 

where 6 and <fi are latitude (zero at the equator; positive northward) and longitude (positive 
in the direction of rotation) in spherical coordinates and in the center-of-mass frame. It 
is found that the angular velocity is insensitive to the latitude and longitude. Thus, only 
the radial dependence in spherical coordinates is considered. If we assume the specific 
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angular momentum, h(m) = |r(m) 2 f2(m), to be conserved during the evolution, the angular 
velocity and surface rotational speed profiles of hydrostatic models can be calculated. 
Figure [JO] shows the specific angular momentum profiles of all considered models in Table EJ 
Since the angular velocity from the FLASH simulations has some variations in the surface 
region, we use the standard deviation of the specific angular momentum in the region of 
0.8 < m/M* < 1.0 to estimate the uncertainty of specific angular momentum. 

Figure [11] shows the evolution of surface rotational speed for all the remnant star 
models. The error bars show the uncertainty of surface rotational speed based on the 
uncertainty from the specific angular momentum. It is found that the rotational speed at 
the surface of all the models significantly decreases to a value less than 30 km/sec. Note 



that t he observed radial velo city of Tycho G is 94 ± 27 km s 1 by 



(120041 ) and 79 ± 2 km s" 1 by 



Kerzendorf et al. 



Ruiz-Lapuente et al. 



(120091 ) . but the linear speed in our models 



are about ~ 200 km s _1 (Table [2]), corresponding to an inclination angle sinz ~ 0.4. Since 
the specific angular momentum h SUTia _ cc oc i? 2 ^/?*) oc R*v TOt , the surface rotation speed 
decreases while the post-impact remnant star is expanding. For the rapidly evolving stars 
such as models A,B, and D, the rotation speed decreases to less than 10 km/sec within 500 
years. After 1,000 — 1,500 yrs, the stars start to contract, slowly increasing the surface 
rotation speed, but the surface rotation speeds of all models at 440 years are still below 
30 km s _1 . In addition, the surface rotation speed of models A, B, C, and D lie below 
10 km/sec at 440 years (Table [3]). Although the surface rot ation speed in our model E is 



still above the upper limit found by 



Kerzendorf et al. 



(120091 ). our progenitor models show 



that the post-impact remnant stars do not need to be fast rotators in the supersoft channel 
(WD+MS). However, Tycho G could also be a stripped giant s tar that has lost most of its 
angular momentum during the SN la impact and then cooled 



an expanded and slowed-down M dwarf as suggested by 



Kerzendorf et al 



Wheeler 



(1201 



I 



20091 ). or 
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Fig. 10. — Specific angular momentum as a function of the mass fraction for all companion 
models in Tabled 
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Fig. 11. — Evolution of surface rotational speed for all the companion models in Table El The 
error is estimated by the variation of specific angular momentum within 0.9 < m/M* < 1.0 
in Figure fTOl 
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CONCLUSIONS 



We have investigated the post-impact evolution of remnant stars in the SDS for SNe la 
via stellar evolution calculations. We examined six possib le binary companion models in 



the mass-orbita! 



Pan et al. 



the setup in 



-period space from 



Hachisu et al 



(l2008bl ) and one companion model in 



(120121). and we pe rformed three-dimensional hydrodynamics simulations using 



Pan et al. 



( 120121 ). The post-impact evolution of the surviving stars was studied 
by using the stellar evolution code MESA, together with the reconstructed hydrostatic 
remnant star models. It is found that the luminosity of post-impact remnant stars increases 
to 10 — 50L Q after the supernova impact and increases to ~ 1OOL within a few thousand 
years, depending on the progenitor model. Due to the energy deposition from the SN 
ejecta, the envelope of the post-impact remnant expands on its local thermal timescale 
(~ 10 2 — 10 4 yrs). After this expansion, stars start to contract and release gravitational 
energy. The post-impact evolution is directly affected by the explosion energy since it is 
related to the amount of unbound mass after the SN impact and the amount and depth 
of energy deposited in the remnant star. Among the calculated models, companion E 
in our simulation (see Tabled]), which has a mass M = 1.44M Q , radius R = A.57R Q , 
effective temperature T e g = 5, 737 K, and luminosity L = 2O.3L after 440 years of 
the SN la explosion, is clo s est to the observed properties of Tycho G as determined by 



Gonzalez Hernandez et al. 



as large as the value given by 



(120091) Although the fits are p r omisin g;, the luminosity is twice 



Gonzalez Hernandez et al 



(120091 ) . Finally, by comparing 



the observed radial velocity to the linear speed in our progenitor models, an inclination 
angle sinz ~ 0.4 can be inferred. The surface rotational speed thus implied (~ 10 — 20 
km s" 1 ) approaches the low upper limit on the rotational speed of Tycho G found by 
Kerzendorf et al.l (120091 ) . Our results provide some support for Tycho G as a possible 
progenitor candidate in the SDS and point to the need for further detailed studies of the 
SDS binary evolutionary channel. 
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